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Abstract 

Within an explicit scalar QED model we compare, at fixed x -C 1, 
the leading-twist K± -dependent 'quark' distribution f q (x, K±) probed 
in deep inelastic scattering and Drell-Yan production, and show that 
the model is consistent with the universality of fq(x,K±). The ex- 
tension of the model from the aligned-jet to the 'symmetric' kinema- 
tical regime reveals interesting properties of the physics of Coulomb 
rescatterings when comparing DIS and DY processes. At small x the 
transverse momentum (kj_) induced by multiple scattering on a sin- 
gle centre is process dependent, as well as the transverse momentum 
broadening occurring in collisions on a finite size nuclear target. 
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1 Introduction 



A significant modification of the quark and gluon distribution functions in 
heavy nuclei - as compared to light targets such as a proton or deuterium - 
is observed in Deep Inelastic Scattering (DIS) experiments on nuclei (for a 
review, see PQI2])- Although such effects manifest themselves on a wide range 
in the light-cone momentum fraction x carried by the parton struck in the 
target, it is useful to discuss separately two relevant kinematical regimes at 
work. In the target rest frame, the typical lifetime for the hadronic fluc- 
tuation of the virtual photon of energy v is given by the coherence length 
l c = 2v/Q 2 = l/(Mx) (where M is the mass of one scattering centre in the 
target, here a nucleon in the nucleus). At large x the coherence length re- 
mains small compared to the typical distance d between two centres, l c < d. 
This is the incoherent regime for which one expects factorization between 
the hard production process on a given centre and the subsequent final state 
interaction of the hadronic (or partonic) fluctuation. Conversely, at small 
x < 1, the fluctuation scatters coherently on several scattering centres. In 
this regime, rescatterings actually affect the hard process. Taking the aver- 
age distance d ~ 2 fm between two nucleons in heavy nuclei, one expects the 
onset of coherence effects such as shadowing below x ~ 0.1. This is indeed 
seen in the small x measurements of nuclear structure functions performed at 
CERN and Fermilab by the NMC and E665 experiments respectively [3]. In 
addition to DIS, evidence for shadowing corrections also comes from Drell- 
Yan (DY) data at small x^ = x measured by the fixed-target experiments 
E772 and E866/NuSea More data are expected at smaller x from the 
Relativistic Heavy Ion Collider (RHIC) facility and in a few years from the 
Large Hadron Collider (LHC). 

Although the nuclear parton densities probed in DIS or through the DY 
mechanism appear to be similar within experimental errors, other (intrin- 
sically non-perturbative) observables depend significantly on the considered 
hard process. In particular, the transverse momentum nuclear broadening of 
the DY pair 01 15] is found to be a factor of 5 or more 1 smaller than what 
is measured in dijet photo- or hadro-production .6]. Note that a similar 

1 Given the various incident reaction energies and depending on the precise definition 
adopted for the nuclear broadening, telling precisely how big the discrepancy is turns out 
to be somewhat delicate although it is statistically significant. 
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discrepancy between DY pair and heavy quarkonium transverse momentum 
broadening has also been reported jH El El • 

In the incoherent regime, one would expect the broadening for large nuclei 
to be proportional to the path length covered by the quark produced in 
the hard process, (pj_) = A 2 A 1 / 3 , coming from the diffusion in transverse 
momentum space due to multiple scattering jS]. Moreover the factorization 
between the hard process and the rescatterings should make the strength 
of the nuclear broadening A a universal quantity (up to trivial geometrical 
and color factors). Yet the smallness of nuclear broadening in the DY process 
seems to contradict universality, which remains so far not fully understood 
El E] • Given the fact that most of the data lie at the borderline between 
the incoherent and coherent regimes, we may wonder whether coherent effects 
could be at the origin of this observation. 

Since rescattering affects the hard process as soon as l c > d, there is a priori 
no reason to expect transverse momentum broadening to be universal at 
small x. However, when comparing two different processes, it is difficult to 
foresee in which process the broadening will be the largest. 

In this paper those questions are addressed by studying in parallel DIS and 
DY within an explicit scalar QED (SQED) model in the small x limit. In 
this limit we will consider two kinematical regimes, namely the aligned jet 
region where the largest component 2u of the incoming light-cone momen- 
tum is mostly transferred to a single final state particle, and the symmetric 
region where it is transferred to two final state particles. In the aligned 
jet region (where soft rescatterings contribute to the total cross section to 
leading- twist), the model is shown to be consistent with the universality of 
the -fT^-dependent distribution f q / T (x, K±) of the target quark participating 
to the hard subprocess, for an arbitrary number of scattering centres in the 
target, i.e., f£jtp(x, K±) = K±). In the symmetric region (where soft 

rescatterings contribute to higher-twist to the cross section), the distribution 
in the Coulomb transfer k± to the outgoing two-particle system is different in 
DIS and DY (for i < 1). This non-universality appears both for a pointlike 
or finite size target. In particular, in the symmetric kinematics the nuclear 
fc^-broadening is driven by monopole rescattering in DY production and by 
dipole rescattering in DIS. 

The paper is organized as follows. In Section |21 we recall the models of 
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Refs. [121 US] for the leading-twist quark distributions probed in DIS and 
DY production on a pointlike heavy target, which we extend to the case of 
a finite size target in Section 12.21 Section |3] is devoted to the extension of 
the model to the symmetric kinematical regime. A summary of our results 
is given in Section HJ 

2 Model for leading-twist DIS and DY quark 
distributions: aligned-jet kinematics 

2.1 Single scattering centre 

Perturbative model for DIS 

The leading-twist i^-dependent quark distributions f q /r(x, K±) in DIS and 
DY are modelled within a scalar QED (SQED) model. Let us start by briefly 
recalling the features of the model for DIS 12j. The contribution to the 
DIS cross section cfdis ( or to the forward DIS amplitude) studied in [T2] is 
obtained by squaring the DIS production amplitude shown in Fig. la. The 
target T is chosen to be a scalar 'heavy quark' of momentum p and mass 
M. The incoming virtual photon of momentum q couples to scalar 'light 
quarks' of mass m, which appear with on-shell momenta p\ and p<i in the 
final state. The 'electromagnetic' charge of the light quarks is denoted by 
e and the light and heavy quarks interact with 'strong' coupling g. We 
work in a target rest frame where 2 q = (—Mxb, q~,0±), xb being fixed and 
q~ = 2u = Q 2 /Mxb — > oo in the Bjorken limit. In this limit the contribution 
to the cross section of transverse virtual photons is subleading (in SQED), 
and we take the photon to be longitudinally polarized in the following. 

Studying the effect of Coulomb soft rescatterings between the light and heavy 
quarks at leading-twist 3 requires concentrating on the aligned-jet region [2], 
where most of the photon energy v is transferred to the struck quark, i.e., 

2 We use the light-cone variables fc ± = fc° ± k z to define a momentum k = (k + , k~ ,k±). 
3 In section 3 we will extend the model to a higher-twist kinematical domain, where 
scales with v, i.e., the ratio y = pT, /{2v) is fixed. 
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Pi ~ q~ ^> P2 ( see Fig- !)■ Moreover, since the hard scale v does not flow in 
the internal propagators of the lower part of Fig. la, the hard vertex 7*q — > q 
is taken at zeroth order in g. The square of the DIS production amplitude 
describes the soft dynamics which can be directly interpreted (apart from 
a trivial factor oc e 2 Q 2 ) as a contribution to the light quark distribution in 
the target f q /T(x, K±). The momentum K corresponds to the quark probed 
in the target and reads K = k — p 2 , where k = ^2 ki is the total Coulomb 
momentum transfer between the light and heavy quarks. In Fig. la we have 
K = pi — q, and the quark distribution is thus probed at x — K + /p + ~ 
—q + /p + = xb- Since K + > 0, it is easy to realize, for instance in a light- 
cone time-ordered formulation, that the hard subprocess in Fig. la is indeed 
7*q — > q, as in the infinite momentum frame. We stress that the diagrams 
contributing to the Q 2 evolution of f q /T being excluded, the model thus 
describes f q /x at an initial soft scale Qq. 

The DIS amplitude of Fig. la has been calculated in ^2] i n the limit xb <C 1. 
In transverse coordinate space, the result including any number of Coulomb 
rescatterings reads 



where we define normalized scattering amplitudes M. in terms of Feynman 
amplitudes M. by 



In Eq. ((H) the factor ijj{r±) denotes the 7* — > qq dipole wave function: 



where a = e 2 /(4ir), y is the light-cone momentum fraction carried away by 
the 'antiquark' p 2 , 




(1) 



*l>(r±) = yVa Q V(m\\ r±) 



(3) 




(4) 



and the function V is given by 





Here r± denotes the dipole size and 



m 



m = p 2 Mxb + m 2 = yQ 2 + m 2 



(6) 
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Figure 1: (a) The DIS production amplitude in the model of Ref. [T2] . 
Coulomb rescatterings are resummed. (b) The DY production amplitude 
obtained from (a) by crossing ^3]. Diagrams where the heavy photon is 
emitted from an internal quark line are suppressed (in covariant gauges) in 
the limit xb < 1. 



The qq dipole scattering amplitude T q q appearing in (JTJ reads 

TWru R±) = ~i ( 1 - e-tf^A-)) , (7) 



where R± is the impact parameter of the outgoing quark and W is the dipole 
single scattering amplitude: 

wfff $ ^_ f d 2 k ± l_^± pin± r kj _ K (\R ± ) - K (X\R ± + rj|) 



(2tt) 2 fci + A 2 - 2 7T 

(8) 

We have introduced a finite photon mass A. Indeed, while the amplitudes W 
and T are infrared safe in the A — > limit 4 , some other quantities are not, 



4 We have for instance liniA^o W(f f ±,R±) = ^ log ( ^^^ J • 
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such as the dipole scattering cross section <r q q(r±) defined below in Eq. (JTUJ). 
Since <J qq {r±) is a basic quantity which will enter our main equations, we 
choose to use the infrared regulator A from now on. This provides a mathe- 
matically well-defined framework, and will moreover allow to compare unam- 
biguously DIS with DY production, for which the infrared sensitivity appears 
at the amplitude level (see Eq. (jTfijl ). 



We stress that 



Tqq(f ± ,Rj 



2lmT q Jr ± ,R A 



(9) 



a unitarity relation which is satisfied by T qq since the latter resums Coulomb 
rescatterings. For later use, we give the scattering cross section of a qq dipole 
of size r±, 



^qq(?±) 









j d 2 R ± 


T m (r±, R±) 





2i j d 2 R ± T qq (f ± ,R ± ) 



(10) 



where the second equality follows from @ and 

d 2 R ± ReT qq (f ± ,R ± ) = 



We give the explicit form of the dipole cross section in our model, 



cr qq( r ±) = 4 / d 2 R±sm 



V 

47T 



K (XR ± )-K (\\R ± + r ± \ 



;i2) 



The DIS production amplitude ((H) incorporates the leading-twist shadowing 
effects discussed in [T2] . and originates from the kinematical domain 

Q 2 , y^oo>P2>^"> h±,Pi±, k~ -> m ^ H 1P2 ~ Mxb ^> Pi ocl/u 

p~ 2 fixed & y = ^ (13) 

In the first line the first and last inequalities arise from the Bjorken limit 
and the aligned-jet kinematics (the latter being stressed in the second line). 
The other inequalities arise from the limit C 1. All scales other than 
v are soft, i.e., intrinsic to the target system. Note that taking a relatively 
large target mass M is not essential to our analysis but will simplify the 
expressions of cross sections. 
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Using dp\jp\ — dp 2 jVz > cross sections will be given by (we will always assume 
V < 1): 

du If d 2 p 2 ± d 2 k± ,-T7,_ r m9 /, ^ 

l^ = ^Jw^^ M(p2± ' k±)l ' (14) 

where the Fourier transform is defined as 

M{p 2 ±M)= [ d 2 r ± d 2 R ± M(r L ,R ± )e- ir ^^- in ^ . (15) 



Model for DY production 



In order to obtain a model for the quark distribution probed in DY produc- 
tion, the first step is simply to exchange the virtual photon and the struck 
quark lines in Fig. la, and to replace q 2 = —Q 2 < by q 2 = Q 2 > 0. The 
momentum of the DY pair has now q + > and reads q = (Mxs,q~ ,q±), 
and the incoming 'antiquark' is chosen with p\± = 0±. We keep the same 
notation for k, p 2 (see Fig. QJi), and K = k — p 2 , implying that K± = q± 
in DY instead of K± = p\± in DIS. In coordinate space the DY production 
amplitude pictured in Fig. lb is simply related to the DIS amplitude by a 
phase factor [T3], 

M DY (r±,R±) = - e ig2G ^ M DIS (f±,R±) (16) 
fd 2 k ± e^_ K (\R L ) 1 / 1 \ 

G ^ = J j^fW^x 2 = ~^r~ ^ ^ 0g {xR-J- (17) 

As already mentioned, in the case of DY production the infrared sensitivity 
shows up at the amplitude level through the phase shift g 2 G(R±). This is a 
direct consequence of the fact that the DY production amplitude of Fig. lb 
involves the scattering of a charge instead of a dipole in DIS. 

Since the phase factor in (|16J1 has no effect on the total DY cross section, the 
present model for DIS and DY is consistent with the universality of the -Re- 
integrated quark distribution. In Ref. [13J it was shown that the non-trivial 
crossing (|T6|) between DIS and DY has nevertheless interesting consequences. 
In particular the distribution da/d 2 k± in the Coulomb transfer k± (k± ^ K±) 
is different in DIS and DY, i.e., non-universal, as we will recall in section 3.2. 

This result is not in contradiction with the universality of the quark distri- 
bution f q /T{x, K±) since in the aligned-jet kinematics (|13|) k±_ is a variable 
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internal to the target structure, integrated out in f q / T (x, K±). Indeed, we 
show now that fq/r(x, K±) is universal within the model of Fig. 1, in agree- 
ment with factorization theorems [TH I15j . The differential DY cross section 
at fixed K± is obtained from (JHJ), 

W iaD ls- = vJlh \Mdy(P^ h = Pk + Ax)| 2 • (18) 
d\ogyd 2 K± 47r J (2tt) z 

Going to transverse coordinate space and using (fTo^) we obtain 

(2tt) 2 daDY _ = ±- f d 2 r ± d 2 R ± d 2 r;d 2 ti;5 2 (r ± + R ± -r^- R[) 
d\ogyd 2 K 1 _ 1» ./ 

Using the constraint from the delta function, we can rewrite the phase dif- 
ference as 

G{R ± )-G{R±) = G(R L )-G(\R ± + r L \)-(G(R^)-G(\R; + r;\)) 
= W(f ± ,R ± )-W(f_[,Rl) , (20) 

since W(f±, R±) = G(R±) - G(\R_l + f±\). It is easy to see that in (P2), the 
phase difference (120)1 can be absorbed in the expression of the DIS amplitude 
given by and ((Zj). After the change of variables r <-> r', R <-> R' we obtain 

( 27r ) 2 dam__ = 1 f d^R^^R^^ + R ± _ ^ _ R^ 

d\ogyd 2 K ± ./ 

x e^- & ^M DIS {?uR^Mms{KM)- (21) 

Interpreting the differential cross sections as f q /r(x, K±) we thus write 

f°Mx, K ± ) = f°/ T s (x, -K ± ) = f°/ T s (x, K ± ) , (22) 

where we used the fact that f^/xix, K±) is a function of K± = \K± \ only. We 
stress that the universality found in ()2*2*|) directly translates into observable 
quantities. Indeed, in DIS K± = pi± (Fig. la) and in DY K± = q± (Fig. lb). 
Thus in the present model the leading-twist DIS struck quark pi± distribution 
and the DY distribution are identical. We stress that in DY the target 
quark distribution is probed at x 2 = K + /p + ~ q + /p + = i e , i.e., x 2 1 in 
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our model. As already mentioned in the Introduction, we expect the coherent 
rescattering physics to come into play at X2 ~ 0.1, where nuclear shadowing 
becomes quantitatively important. 

One might ask whether the universality (J22|) would be preserved in a more re- 
alistic model for DY production with a composite projectile. The role of spec- 
tators in our model for DY production is studied in Appendix A. We find that 
spectator rescatterings do not affect the DY distribution do~DY/d 2 q±_ (see the 
right hand side of Eq. (jHTjl which only depends on the DIS amplitude). The 
main consequence of using a composite projectile is to replace the infrared 
cut-off A (the finite photon mass) by the inverse size S of the projectile. We 
also show that the model with a composite projectile is consistent with fac- 
torization (see and Ref. [T5]). Since spectator rescatterings do not lead 
to any breaking of universality between the DIS and DY ^^-distributions 
(at least in our SQED model), we will neglect (projectile) spectators in the 
following, and use for DY production the model with a pointlike scalar pro- 
jectile presented above. 

It is worth recalling that the K± distribution 5 can be expressed in terms 
of the qq dipole cross section (see for instance [E|). From (|2Tj) and ((TJ) we 
obtain 

(27T) 2 = [ d 2 f ± d 2 f;e-^- r ^- R ^(r ± )^(r;) 

d log y d 2 K± J 

x J d 2 R L T qq (f ± , R ± )T*^,R± + r± - ?l) . (23) 

With the identity 

T q q(fj_, R±)T*^(f^, R± + f±- r_i) = 

iT m (f_L, R±) - iT*^(fl,R x +f±-fl)- iT q q(fj_ - f_[, R±) (24) 
the expression (|2l^l becomes (using (|TU|) ) 

( 27r )2 = f d 2 r L cZ 2 f J _' e - i ^-^')-^V'K)^(^) 



d log y d 2 K 1 

X [|<7qq(r±) + Kq(r±) - ^(\f± - f[ \ )] . (25) 

The latter expression can be found in |16J (in the more general case of finite 

vY 

From now on we suppress the subscript 'DIS' or 'DY' for this universal distribution. 
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2.2 Nuclear target 



The model for DIS and DY production on a pointlike target described in the 
previous section can be directly generalized to the case of several scattering 
centres. Since the DIS and DY amplitudes off a single centre have been 
derived in the limit <C 1, they should also describe the production off a 
'nuclear target' in the limit where the nuclear radius R A is kept smaller than 
the coherence length l c , 

R A «l^ J L-. (26) 

In this limit the production amplitudes are simply obtained from (0) and (|16|) 
by replacing the scattering potential on a single centre G(R±) given in (JT7J) 
by the scattering potential on A centres located at transverse 6 positions X{±, 

G(R±) = J ^ °a(R±) = E g (l & - • ( 27 ) 

Denoting 

W A (f ± ,R ± ) = G A (R ± )-G A (R ± + f ± ) (28) 

T*(?±,R±) = -i (l - e-tfw^fiJ} , (29) 

we obtain for the DIS and DY amplitudes on the nuclear target A 

M^ IS (r±,R±) = s/^^(r ± )T^(fs_,R ± ), (30) 
•M£ y (riA) = -e^ G ^ M^ IS (r ± ,R ± ). (31) 



We can repeat the steps leading from ([T%]) to (J22J to realize that M-dis an d 
TW^y provide a model for the 'nuclear' quark distribution function which is 
consistent with the universality of the -fT^-dependent distribution 

f^(x,K ± ) = f° A s (x,-K ± ) . (32) 

The physical content of f^ A s = in terms of Coulomb rescatterings on 
the A static centres is depicted in Fig. El 

6 In the total coherence limit (|26|l the longitudinal positions of the centres are irrelevant 
and only the thickness function T of the target enters (see below). 
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Figure 2: Model for the DIS (a) and DY (b) nuclear quark distribution 
functions f^/^(x, K±) and f^(x, K±). In DIS we have K± = k±—p 2 ± = pi± 

and in the DY case K± = k±_ — p 2 ± = q±. The model takes into account 
any number of Coulomb rescatterings on every centre. In the total coherence 
limit (|26j) the A centres only differ by their relative transverse positions. 



We now evaluate the universal distribution da A / d 2 K± by averaging over the 
positions of the scattering centres. For our purpose it is sufficient to average 
with a uniform distribution, 



A ' J2; 



< ^imMr) • (33) 



i=l 



where S — J d 2 x is the transverse area of the target. Inserting (jHUjl in (|2*T|) 
we find 

(27T) 2 d ° A _ = f d 2 f ± d 2 rle~ i( - r ^>^tl;(r ± )i;(r;) [ d 2 R ± 



d log y d 2 K± 

x (<^(T ? ± ) i2±)-x7^(f J :,i2 J :)*-i^(f J .-r 1 :.7?x)> i . 



where R[ = R± + f± - r[. Using W A (r±,R±) = Ylt=iW(r±,R± - x iA _) 
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and the approximation 7 <J q q{ r ±) S (valid when the nuclear radius is much 
larger than the interaction range A -1 ), we show using (|1U|) that 

l-iT A (f ± ,R ± )) =0 (35) 



' qq 

for R± outside the nucleus, and 

for R± inside the nucleus, where T = A/S is the target 'thickness'. The 
average (j3Hj) of the qq dipole matrix element S q q = 1 — iT qq - resulting in the 
form (|36|) allows to interpret the rescattering process as a random walk in 
the transverse plane (see Ref. [TT]). 



r j n A 

(2tt) 2 — — = S I <Pr ± d l rle- i( ? ± - p ±> K - L i/>(r ± ) i[>(r A 



The equation (J3*3|) becomes 
, da A 
dlogy d 2 K 

I _ e -T^(r ± )/2 _ e -T CTq<5 (r ± ')/2 + e -Ta m {\f ± -r [ |)/2 1 _ ^J) 

This corresponds to an 'eikonalisation' of the distribution off a single scat- 
tering centre (J25|) (see also ^Ej). In the limit of vanishing thickness T — > 0, 
(jnij) reproduces (j25j) up to the factor A. Since da A /d 2 K± is universal, we see 
that the DY ^-distribution (and a fortiori the total DY cross section) can 
be expressed in terms of the qq dipole scattering cross section, as is the case 
for DIS OSj. 



3 Extension of the model to the symmetric 
kinematical regime 

3.1 Symmetric kinematics and interpretation 

In the preceding section we have focused on leading-twist contributions to the 
DIS and DY cross sections, arising from the aligned-jet region defined in (|13|). 

7 With this approximation and the uniform distribution l|33|) , we can derive Glauber- like 
expressions without resorting to the full formalism of [17| . 
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We would like to stress however that the validity of our calculations so far 
is not restricted to the aligned-jet domain. Since we only used p 2 ^ u i 
they remain correct even when p 2 scales with v (i.e., in the 'symmetric' 
kinematical region), provided the ratio y = p^/(2u) is fixed to a small finite 
value, i.e., in the region 

p~ ~ v — > oo > p 2 > M > ki±,p i± , fcr, m > ~ Mi B > , oc l/u 

fixed y = — -C 1 -v^ £>2 scales as z/ . (38) 

In fact our results can easily be extended to the domain where y is not small 
as compared to 1. For simplicity we concentrate on the n « 1 limit, but 
expect our following considerations to apply also to situations where y ~ 1, 
such as quarkonium (y ~ 1/2) and dijet leptoproduction. 

With the symmetric kinematics (|3*Hj) the interpretation of the contributions 
to the DIS and DY cross sections we evaluated is modified. 

(i) The antiquark of final momentum p 2 is now part of the hard subprocess, 
which reads 7*g — > qq in DIS (Fig. la) and qg — > qy* in DY (Fig. lb). 
The DIS and DY processes are now interpreted respectively as dijet 
leptoproduction and the associated production of a DY pair and a jet. 
In the symmetric kinematics these processes probe the /cj_-dependent 
'gluon' distribution / g /r(^B,fc±)- 

(ii) The transverse momentum transfer to the hard system now corresponds 
to the Coulomb k± exchange, i.e., not to K± = k± —f>2± any longer. In 
the case of DIS, k± = pi± + p 2 ± is the momentum imbalance between 
the jets. In DY, k± — q± + f>2± is the imbalance between the DY pair 
and the associated produced jet. 

In the following we concentrate, in the region (}3*8*j) . on the DIS and DY k± 
distributions, and show that those are different. 

In section 3.2 we show that the fc^-distribution on a single scattering centre 
becomes non-universal beyond leading order in g 2 . In particular higher order 
corrections to daoY I 'd 2 k± vanish in the A — ► limit (at fixed k±), contrary to 
DIS. In section 3.3 we use the model for a nuclear target presented in section 
2.2 to investigate the target size dependence of da A /d 2 k±, and find that the 
strength of fcx-broadening depends on the process. 
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3.2 Production off a single centre: non-universality of 
the /^-dependent distribution 

In this section we show that in the symmetric kinematics (|38J) the DIS (dijet) 
and DY (+ jet) /^-distributions are different. We first consider the DIS 
model presented in section 2.1. From (|14|) the differential DIS cross section 
at fixed k± reads 




(39) 



where 




(40) 



From we obtain 




(41) 



where we defined 



da q <i(r ± , kj_) _ |T q q(rl, k ± )\ 2 
d 2k ± ~ (2tt) 2 



(42) 



Let us turn to the case of the DY process, 




(43) 



where -M_Dy(rj_, kj_) can be shown using (JTJ) and (fTH|) to be proportional to 
its own Born value, 




(44) 



C(k 2 ± ) M^T(r±,k ± ) , 



(45) 
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where we define 



C(k] 



T*° m (?±, k, 



fcj + A 2 
ig 2 



d 2 R l e ^ 2 G(R^) e -iR^ 



-V^(r ± )T^(f ± ,k 



0^2 sin (^- fc i-/ 2 ) 



A; 2 +A 2 



(46) 
(47) 
(48) 



Using (jUJ), and (|42j) we find 



dlogy d 2 k± 
This can be reexpressed as 



\C(ki)\' d'f ± mr ± )\ 



d?k\ 



(49) 



(27T) 



da 



d log yd 2 A; j 
by using (|47)|) and the identity 

Using and (joTi|) we also obtain 

2 \|2 



d 2 rHVK)| 2 4sin 2 (^-^ w ' 2r " VA r 1 



(50) 



d 2 6e^ / d 2 /? ± (*T q * q (6, i? ± ) + 1) . (51) 



\C(ki)\ 



Jd 2 be^a m (b) 
Jd 2 be^a^ Tn (b) 



(52) 



Comparing (|4*T|) and ()49|) it is clear that the fci-distribution is process depen- 
dent. We can stress this point by noting that for k±^> X, the DY distribution 
(|4T?|) equals the Born distribution, contrary to the DIS distribution. Let us 
consider the kj_ — > oo limit at fixed A. In the expression of IC^Ar^)) 2 given in 
(|52|) . the phase factor e lb ' k± rapidly oscillates, except if b < l/kj_ — ► 0. When 
6 — > the dipole cross section can be approximated by its Born value (see 
(JT2J)). Thus when k± > A we have \C(k 2 ± )\ 2 ~ 1, yielding 



da 



DY 



d log yd 2 k j 



fc±>A 



d 2 f± \ip(r±] 



da*°™(f ± ,k A 
d 2 k> 



da%f n 
d log yd 2 k j 



(53) 
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We can obtain this latter equation more rigorously by using (]4*5]l and cal- 
culating C(kj_) given in (|46|) for A — > 0, the other scales being fixed. Using 
G(R±) ~ -log(A J R ± )/(2vr) we get 

C(kl) xZo -^^e^( log (^)^) [1 + (A)] =► \C(kl)\ = 1. (54) 

Thus C(A; 2 A 2 ) is a pure phase factor, leading to (]53Jl . which confirms to 
all orders the result shown in at next-to- leading order in g 2 . 

Since the k±-integrated DIS and DY cross sections are identical (and different 
from their Born value), the result ([53]) implies that in the A — > limit 
and in the DY case, only vanishing k± ~ A — > contributes to Act = 
a tot _ cr Bom ^his difference observed between the fc^-dependent DIS and DY 
cross sections is similar to what Bethe and Maximon found in the case of high 
energy pair production and bremsstrahlung [15] . In the present context the 
effect clearly arises from the infrared divergent Coulomb phase (for A — > 0) 
in the DY production amplitude, and suggests that in collisions on a proton, 
the k± exchange in DY + jet production might be smaller than in dijet 
leptoproduction, and thus not 'intrinsic' to the target. 

In a realistic situation, when the projectile and target are composite and 
neutral, the DY production amplitude is infrared finite. Effectively, the role 
of the infrared cut-off A is played by the largest infrared momentum cut-off at 
disposal, given by the inverse size 8 of the smallest incoming hadron, projec- 
tile or target. We expect the typical transverse momentum contributing to 
Actdy to be k\ ~ 8 2 instead of k\ ~ A 2 , as discussed in the end of Appendix 
A. 

In the next section, we investigate the non-universality of the /c^-distribution 
in the case of a nuclear target. 

3.3 Nuclear target: non- universality of /c^-broadening 

The 'nuclear target' model of section 12.21 is now used to derive the target 
size dependence of the DIS and DY /^-distributions, where k± is the total 
transverse Coulomb exchange. We have shown that the distribution in K± = 
k±—p2± is universal (see ([32]) ). and this holds in the two kinematical domains 
(H3J and 
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On the contrary, we now explicitly show that the distribution in the trans- 
verse Coulomb exchange da A /d 2 k± is process dependent. Similarly to the 
case of a single scattering centre, this can easily be guessed from Eq. (p?T|). 
The DIS and DY amplitudes differ by a pure phase factor in R± space, but 
there is no reason to expect this to hold in the space of the conjugate variable 
k ± . ' 

Inserting (gUj) in (gHj) we get 

d °Dis r^u, m2 j^g^jjj 

- ' d r ± \ip{r±)\ -= , (55) 



d\ogyd 2 k± J d?k 

with 

da A = JL / d 2 R ± d 2 R; e -^- R ^ x (56) 



A 



d 2 k ± (2vr 

^1 _ e -ig 2 W A (f ± ,R ± ) _ e ig 2 W A {r^,R[) + e -ig 2 {W A (r^,Ri_)-W A {f^Rl))\^ 



representing the differential elastic scattering cross section of a qq dipole of 
size r± on a nuclear target. The form (JHH)) emphasizes the decoupling between 
production and rescattering of the qq dipole in the process. The average over 
the positions of the scattering centres leads to 



da± S 



d 2 k ± (27T) 2 



Using we eventually obtain 

dn A 
d\ogyd 2 k 



(27T) 2 da ° Is ^ = S I d 2 r ± \iP(r ± )\ 2 I d 2 be^ 



jl - 2e-^x)/2 + e -T/^m-e^)*^j (5g) 

A similar calculation for DY production is performed inserting the Fourier 
transform of (J3TJ) into (|43|). 

(27T) 2 ^ 



d log y d 2 k j 
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giving 



d 2 r ± \ij(r ± )\ 2 
d 2 rMri)\ 2 



dCTpy 



d 2 R\ ( e i92GA(R ^ - e^GAd^+r^n e ~iR±-k 



'l-e l 



' d 2ft ig 2 G A (R±)-iR±-kj 



d<?t 



dlogy d 2 kj_ d 2 k\ 



d 2 fj^(r ± )| 2 4sin 2 (^-^ 



.4 

(59) 
(60) 



where da^/d 2 k± is the differential elastic scattering cross section of a quark 
on a nuclear target in our model, i.e. 8 : 



d<?q 



d 2 k l 



1 



fcx^O 



(27T) 



d 2 R<d 2 R:e-^- n ^ ( e -if(G A & ± )-G A (R±)) 



(61) 

We perform the average using Ga(R'±) — Ga(R±) = Wa(R± — R±, R'±) and 



dai 



d 2 ki 



S 



d 2 be^ e~ T ^ b)/2 



(62) 



finally leading to 
da™ 



(2 



7T 



'DY 



dlogy d 2 k\ 



S I d 2 r ± \*lj(r ± )\Hsm 2 



d 2 be^^ e ~ Tc7 ^ (b)/2 . 

(63) 



In line with the above discussion for DIS, exhibits some decoupling 
between production and rescattering, up to the factor 4 sin 2 (f± ■ k±/2) spe- 
cific to DY production. Notice that despite the formal appearance of cx q q(6) 
in its expression ([62)1. da^/d 2 k± represents monopole elastic scattering 9 , in 
contradistinction with the dipole scattering cross section (JSTj) appearing in 
DIS. 

Comparing (J55)l and (|50|l we interpret the non-universality of da A /d 2 k± in 
our model as a direct consequence of the type of object which interacts with 

8 We neglect contributions ~ S^ 2S> (k±) which do not contribute to da^ Y /d 2 k± due to 
the factor s'm 2 (f ± ■ k±/2) in i|5U|) . 
9 This point is discussed in detail in 



18 



the nuclear target. In the rj_-integral giving d<r A , IS /d 2 k^, the dipole wave 
function if) selects r± < l/m\\ in da A ^/d 2 k±. Therefore, the hard scale 10 
mil enters the physics of rescattering and is expected to play a major role 
for k± broadening. Conversely, in (p)Uj) the scale m\\ enters d<T A , Y /d 2 k± only 
through a target-independent factor and is therefore not expected to govern 
k± broadening. 

In order to display the differences between DIS and DY we investigate the 
ratio 

where a p and a A denote the DIS or DY cross sections off a single scattering 
centre and on A centres, to leading order in g 2 and next-to-leading order in 
the target thickness 11 T = A/S. By expanding (|58|l and we get, for 
A <C kj_ <C m\\ and keeping only the leading logarithms: 

R D y{k±) = 1 + l0g^ + O(T 2 ). (66) 



The latter results explicitly demonstrate the non-universality of the /in- 
dependent distributions which was already apparent when comparing the 
full expressions (J55|) and ()6(J|). As compared to the production on a single 
centre, and in the region under consideration A <C k± -C Tom, dcr A /d 2 k± is 
slightly reduced in DIS and strongly enhanced in DY production. In the 
DY case, we recall that for a single scattering centre (see section 3.2), the 
k± distribution for k± ^ X equals that at Born level, i.e., in the absence of 
rescattering (see Eq. (|55|l). This might have suggested /^-broadening to be 
reduced in DY as compared to DIS. The explicit calculation with a finite size 
target of 'thickness' T shows that this does not happen: the deviation from 
unity of the ratio -Rdy(&±) survives the A — > limit. 

10 In the symmetric kinematics (13811 we have m?i = yQ 2 + m 2 ~ yQ 2 . 

n Note that the limit of small target thickness T oc A 1 / 3, is consistent with the total 
coherence limit Ra <C l c (see (J3BJ), in which our DIS and DY amplitudes l|50|l and (|3T|l 
have been derived. 
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It is simple to realize that the fc^-broadening defined as 

A(kl) = (ki) A -(kl) p , (67) 

although different in DIS and DY, scales as g 4 T both in DIS and DY. Evalua- 
ting A(k\) from and (and using do~ p /d 2 k± oc l/k±) puts some light 
on the difference between DIS and DY. In DIS the small probability ~ T/m| 
for the rescattering of a dipole of size ~ l/ m || * s compensated by a large 
typical momentum transfer k\ ~ rajj. In DY the hard scale my does not enter 
the expression of R(k±), and the relatively large (monopole) rescattering 
probability ~ T/k\ is now compensated by a small typical transfer k\ -C mij. 
As already discussed, the non-universality of nuclear fc^-broadening in the 
coherent limit studied in our model is a natural consequence of the type of 
object (dipole or monopole) which interacts with the target. 



4 Summary and outlook 

Within an explicit scalar QED model, we have studied transverse momentum 
distributions in the coherent limit x <C 1 for DIS and the DY process. In the 
aligned jet kinematics, where the leading quark or the DY pair carries most 
of the projectile momentum, the distribution in the transverse momentum 
K± of these particles is universal, both for pointlike and extended targets. 
This is consistent with the universality of the K±- dependent target quark 
distribution. On the contrary, in the symmetric kinematical region, the rel- 
evant transverse momentum k± (k± ^ K±) is that of the hard subsystem, 
i.e., of the quark-antiquark pair of the photon fluctuation in the case of DIS, 
and of the DY virtual photon and the final antiquark in the case of DY. The 
transverse momentum transfer k± between the target and the hard subsys- 
tem is different in DIS and DY already for pointlike targets. The extension to 
a finite size 'nuclear' target stresses the physical origin of this difference. In 
DIS the qq dipole rescatters with a small probability but undergoes large k±- 
kicks, whereas the DY k±_ distribution is sensitive to monopole rescattering, 
more likely but involving smaller kicks. 

We stress that our nuclear transverse momentum distributions are expressed 
in terms of the dipole cross section cr q q and of the thickness function T, and 
contain factors of the type 1— e~ Tcr which are expected for classical scattering. 
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This behaviour is a consequence of the statistical average we have performed 
on the positions of the scattering centres. Thus also in the coherent region 
rescattering off an extended target turns out to have the nature of stochastic 
multiple scattering [20] ■ Coherence is important in that it fixes the nature 
of the object (dipole or monopole) that rescatters. 

The non-universality of the k± distribution found in our scalar QED model 
in the coherent small x limit calls for a systematic study of the nuclear k±- 
broadening measured presently in fixed-target experiments (x ~ 0.1) or in 
the near future at RHIC (x ~ 0.01). Indeed, at the present stage we cannot 
give a quantitative answer to the puzzling observation mentioned in the In- 
troduction, namely the smallness of transverse momentum broadening in DY 
production, as compared for instance to the broadening of the dijet momen- 
tum imbalance in dijet photoproduction. We however emphasize that with 
our notations the observed transverse momentum is q± = K± in DY produc- 
tion instead of k±_ in dijet photoproduction (analogous to our DIS process in 
the symmetric kinematics). We have shown that daoisi 'd 2 k± ^ daDY/d 2 k±, 
thus it is even more natural to expect the distributions daois/d 2 k± and 
doDY I 'd 2 q± with respect to distinct transverse momentum variables to be 
different in the coherent regime. This can be actually checked explicitly in 
our SQED model. The universal ratio 



can be obtained from (|37j) and 1)25)1 . and reads at leading order in g 2 and for 
K ± < m\\\ 



We see that contrary to the 'dijet' leptoproduction ratio B,Dis{k±) given in 
(jHSJ), the ratio Rdy(q± ^ m \\) exceeds unity. Small g_i_'s are favoured in DY 
production off a nucleus. Whether the latter result, obtained in our abelian 
model, can explain the observed smallness of g^-broadening in the hadronic 
world will be addressed in a future work. 
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(68) 



2a 1 /mil \ 

R(K ± ) = R DY (q ± ) = R DIS (Pi±) = 1 + ^Hr lo § 2 ( -r ) + ( T ') • ( 69 ) 

unit \ A / 
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Appendix 



A DY production with composite projectile 

Here we study the role of spectators in DY production. For this purpose 
we extend the model of Fig. lb to the case of a composite projectile. The 
corresponding model is depicted in Fig. El 

The fluctuation of the 'hadron' projectile of mass and momentum into 
the spectator quark (of final momentum p Q ) and the active antiquark (of final 
momentum p 2 ), is described by a scalar cubic coupling ig . The projectile is 
chosen to carry the 'electromagnetic' charge e but to be neutral with respect 
to the 'strong' interaction of coupling g. The spectator quark has no electric 
charge but strong charge g. The large incoming light-cone momentum P^ 
splits into the active antiquark and spectator momenta with finite fractions 
z and 1 — z. The different momenta appearing in Fig. 3 read (we choose 
mh < 2m in order to forbid the h — > qq decay): 

2 

P h = (-JL,2u 7 ± ) ; p=(M,M,0 ± ) 

,pl ± + m 2 _ sPh+m 2 _ _ . 

P° = (n vT' V 1 _ z ) 2v >Pa±) 5 Pa = ( = >P2>P2±) 

(1 - z)2u p 2 

q = ( Q z2 ^ Q±) - ( Mx b, z2u, q ± ) , (70) 




Figure 3: Model for DY production with a composite projectile. The spec- 
tator is produced with final momentum p Q . 
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where now xb = Q 2 /(z2Mu). We will use again the limit <C 1, as well as 
the kinematics defined in 

The DY virtual photon can be radiated either by the projectile or by the 
active antiquark, both having electric charge e. However, in the kinematics 
(Unj), the typical times associated to the fluctuations h — > qq and q — > 7*q are 
respectively of order vjp\ L — > oo and 1/Mxb- Thus the photon is radiated 
after the h — > qq fluctuation. The diagrams where the virtual photon is 
emitted from the projectile are suppressed (in Feynman gauge) by a factor 
~ pl ± /(isMx B ) ~ O (pz jv) according to (fT3|l . 



A.l Consistency with factorization 

The covariant calculation of the DY production amplitude of Fig.|3]is similar 
to the calculations performed in Refs. [T2*l ll3|. The leading-twist contribution 
is obtained for a virtual photon with longitudinal polarization 

=(-,--, ± ). (71) 

zv zv 

Going to transverse coordinate space, 

M{r u R ± , u ± ) = J M(p 2± , k ± ,p 0± ) ^-^fti) 

(72) 

and resumming Coulomb scatterings yields the result 

M D y(r ± , R ±t u x ) = M DY (r ± , R ± ) e -^\R M) 0(^0 (?3) 

z 

Here M.dy(^±, R±) is the DY production amplitude in the absence of specta- 
tor obtained from (J2J) and (fTSjl. the function G is defined in (fTTj) and <f>(z, u±) 
is the h — > qq wave function 

= <? z(l (74) 
5 2 = m 2 _ z (i _ z )ml (75) 

which can be represented as 

<K*,u±) = /|fe^Po±)e^^ (76) 
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z(l — z) 

<P{z,Po±) = 9o p2 - p ■ (77) 

The phase factor in ()73|) arises from Coulomb rescatterings of the spectator 
quark. Those indeed contribute to the DY production amplitude of Fig. El 
since vanishingly small light-cone energies kf oc 1/v can be transferred to 
the spectator without any cost. The finite energy k + — ~ O(Mxb) 

is transferred to the active antiquark in order to produce the final state 
invariant mass ~ Q 2 . The phase in ()73|) is infrared divergent (the spectator 
carries the charge g) but this divergence compensates that appearing in (fTtj)). 
as expected for dipole rescattering 



M DY (f ± , R ± , u ± ) = - e ^C^) M DJS (f ± , R±) , (78) 

z 



where we used G(Rj_) - G(\R± +u±\) = W(u±, R±). 

We now proceed as in section 2 (see Eq. (fTHj) and following). In the presence 
of the spectator the differential DY cross section is of the form 12 

da DY f d 2 p 2 ± d 2 p ± 



d 2 qi 



/d f>2± d po± -, -, mo 

T^TTp" "(2vrj2 ' M dy{P2±, k ± = p 2 ±+Po± + q±,Po±)\ . (79) 



Going to transverse coordinate space and using (|73*|l leads to 
da by 



d 2 K\ 



oc j d 2 r±d 2 R±d 2 u±d 2 r^d 2 R' ± d 2 u' ± (fi(z,u±)(f)(^ 

x 5 2 (f± + R±-r±- R_i_) 5 2 (u ± + R±-u_i~ R±) 

x e ig 2 (G(R ± )-G(R^)) e -i(G(\R ± +u ± \)-G(\Rl+ul\)) 

x e~^- n ^ M ms (r ± , R±)M* DIS {rl, R±) . (80) 

From the (^-constraints the Coulomb phase associated to spectator rescatter- 
ing cancels out, G(\R± + u±\) — G(\R ± + u±\) — > 0, and the remaining phase 
difference G(R±) — G(R ± ) is absorbed in the expression of AAdis given by 
(JTJ) and ((Zj), leading to 

^ oc / d 2 r L d 2 R L d 2 rid 2 Rl 6 2 (f ± + R ± - r[ - R[) 
d 2 q± J 



12 For the purposes of the present Appendix we do not need to specify the normalization 
of differential cross sections in the following. 
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x e-^-flxKx M DIS (r x ,Ri_)M* DIS (K,Rl) 

x J d 2 u_ L 4>(z,u±)4>(z,\u ± + R_l- R±\)* . (81) 

Thus spectator Coulomb rescattering does not affect do by I 'd 2 q± (and a for- 
tiori not the total leading-twist DY cross section either), which is consistent 
with factorization, as we briefly see now. 

From (fTBj) one gets 

J d 2 u^(z,u ± )<P(z,\u± + R±-R:\y = I ^l^poi^e-^-*^ 

(82) 

Inserting this into (|%T|) and identifying J d 2 p ±\ip(z , Po±)\ 2 with the projectile 
antiquark distribution fq/h(z, —Pol) we obtain 



dopy 
d 2 q ± 



oc J d 2 p 0± fq /h (z,-po±)f q / T {xB,Po± + q±)- (83) 



The latter equation shows that our DY model with spectator is consistent 
with factorization theorems involving i^^-dependent parton distributions 



A. 2 Non-universality of k± Coulomb exchange 

Here we argue that the result found in section 3.2, namely that the typical 
k± contributing to Aody = o'dy ~ ^By 11 is /c± ~ A ^ 0, naturally translates 
to k± ~ 5 in the case of a composite projectile of size Rh ~ 1/5. 

The fc^-distribution reads 
da DY f d 2 p 2 ± d 2 p Q± 



d 2 k ± 

oc 



f d 2 p 2 ± dy ± ... , w ? - x,2 



x gig (W(u±,Rx)—W(u±,R ± )) M DIS (f±, R±)M DIS (f ± , Rl) , 

(84) 
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where we used ()78|) . The integrand of (jHjj) depends on the scales m\\ and 
S, the latter corresponding (see (173)1) to the inverse transverse size of the 
projectile. When 5 — > 0, u± ~ 1/5 —>■ oo, and the (finite) Coulomb phase 
in (1S41 becomes 



W(u ± , R ± ) - W{u ±J R_i) = G(R ± ) - G(R_i) + log ' 1 " L B 



2tt \|ul + .R/| 
- G(iO (85) 



We thus recover, in the 5^0 limit, the fcj_-distribution (|43j) in the DY 
model without spectator, for which we have shown that k± ~ A — ► 0. In 
other words, when the size of the projectile R^ ~ 1/5 — > oo, the spectator 
plays no screening role any longer. In practice 5 is non-zero, 5 ~ m (but 
still 5 <C my ~ \fyQ in the kinematical region (jSHJ)) and the typical 
contributing to Aa^y is of order 5, the largest infrared cut-off at disposal. 
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